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In statistical physics the mean-field approximation is 
one of the most common and easy-to-use frameworks. It 
is also one of the most powerful, and often the only avail- 
able one. It allows one to convert the study of a many- 
body problem of interacting degrees of freedom into that 
of independent degrees of freedom. There are several 
manners of performing a mean-field approximation. Our 
purpose is not to review them nor to discuss the cor- 
responding well-documented pitfalls 0, 0, 0, 0> 0> El- 
Among those various mean-field versions, the one that 
is based on a rigorous variational principle plays a spe- 
cial role. Our purpose in the present note is to confront 
this rigorous, albeit cumbersome approach, with simpler 
and widely used formulations. We would like to analyze 
a hazardous ambiguity in the concept of free energy in 
those apparently more physical formulations which, to 
the best of our knowledge, has not been noted before. 
The reason why these conceptually erroneous statements 
that can be found in many places in the literature (@|, 
) have never been challenged is that pedagogical pre- 
sentations are usually confined to the ferromagnetic Ising 
model, which has a scalar order parameter, and for which, 
somewhat luckily, the dangers we will point out remain 
hidden. 

We have chosen to illustrate our discussion with the 
anti-ferromagnetic Ising model on a two-dimensional 
square lattice of N sites, say with periodic boundary con- 
ditions for our discussion would be pointless with a 
scalar order parameter. The N spins s, = ±l,i = l,...,iV 
are interacting according to the following Hamiltonian: 



H = 



J ^ SiSj 



(1) 



where the sum Yluj) runs over the 2N distinct pairs of 
nearest neighbour sites, and where the constant J > 
is the anti-ferromagnetic coupling. The high temper- 
ature phase of the system is paramagnetic (rrij = (si) 
vanishes). As the temperature is decreased below the 
so called Neel temperature Tjv anti-ferromagnetic order 
sets in: The spins align in opposite directions on two 



square sub-lattices, as shown in Figure ^ The lattice 
sites are conveniently divided into two subsets a and 6, 
as already depicted on Figure ^ 

We begin with a reminder of the variational formula- 
tion of the mean-field approximation (route 1), and then 
present a more standard approach (route 2) and the ac- 
companying difficulties in interpreting the related free 
energy. 

Route 1: Variational procedure. The idea |9l llfl ITT| is 
to introduce a trial Hamiltonian Tt m i depending on the 
original degrees of freedom {si} and on two parameters 
m a and mb the physical meaning of which will become 
clear later. An intuitive choice for 7i m f is 



Hmf = +4Jm h ^ Si + 4 Jm a ^ , 



Then one splits TL into 



W = H m t + (H — H m {) 
so that the free energy F of the system reads 
F = F mf -fc B Tln(e-««-*->} mf . 



(2) 



(3) 



(4) 



Here, /3F m f — — logZ m f where Z m { is the partition func- 
tion associated to (3 — l/(fc B T) is the inverse temper- 
ature with k B the Boltzmann constant and the angular 
brackets (...} m f denote an average using the Gibbs mea- 
sure related to 7i m f, that is with weight exp(— (37i m t). 
The exponential being convex, one is led to the inequal- 
ity 



F < F mi +(H- H n 



(5) 



Equation JSJ) often appears under the name of Bogoli- 
ubov inequality and may be used to find the best set 
of parameters m a and mt that render (f>(m a ,mb) — 
F m {+('H—'H m {) m { minimum, that is as close as possible to 
the exact free energy F 0| . It is important to note here 
that the best approximation for F is not F m {(m a ,mb) 
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but <p(m a ,mb). Using that (sigo)mf — — tanh4/?Jmb and 
( s jeb)mi = — tanh 4/3 J?n a we arrive at 

/3F mt 



jV r 

(3(p(m a ,mb) — — — In (4 cosh 4/3 Jm Q cosh 4/3 Jm?,) 
+ 4/3 J tanh 4/3 Jm a tanh 4/3 Jm fc 
+ 4/3 J (m a tanh 4/3 Jm a + mt tanh 4/3 Jm;,) 

(6) 

Extremizing 6 leads to the set of equations 



86 

= 0, ^- = 



(7) 



9m a 9m{, 
==> m a = — tanh4/3Jmfc, = — tanh4/3Jm Q (8) 

The latter system of equations has a unique solution 
m a = mi, = at (3 < /3jv = jj (TV = 4J/fc B ) and 
possesses an additional set of two nonzero solutions for 
/3 > Pn {Tn is identified as the Neel temperature). In 
the high temperature phase the paramagnetic solution 
m a = mt = becomes the global minimum of 6, just 
as the nonzero solution m a = —nib =/= does in the low 
temperature phase (one can verify that the matrix of the 
second derivatives of 6 is positive definite at those ex- 
trema). Note also that for (3 > (3 m the paramagnetic 
state is a saddle point of 6 with the unstable direction 
along the m a — —nib line (see Figure^. Right at the 
minimum, the expression of 6 reads 



6(m a ,m b ) = Fmf- 2NJm a m b 



(9) 



where m a and mb are the solutions to the system in 10. 

At this stage we have simply postulated a trial Hamil- 
tonian H m { without providing much of a physical motiva- 
tion. It is a posteriori clear that H m i describes a system 
of independent spins in an external magnetic field. For 
spin Si of sub-lattice a this magnetic field is interpreted 
as the mean magnetization resulting from the four near- 
est neighbors on sub-lattice b, as is confirmed by the fact 
that at the minimum of 6 one can indeed verify that 



\Sea/mf : 



tanh4/3Jm b = m a , {s eb ) m { = m b (10) 



In practice, however the variational procedure is not 
physically transparent and is mathematically rather 
heavy. Furthermore it must be supplemented with a 
reasonable input of physical intuition when postulating 
a trial Hamiltonian, lest the outcome of the calculation 
should be dull. Hence, for all these reasons, in spite of 
6(m a ,mb) being a bona fide mean-field free energy, it is 
rarely used in standard courses. The purpose of the se- 
quel is to present an alternative and widely used formula- 
tion of the mean-field approximation 0,0) which at first 
glance appears more satisfactory on physical grounds, 
but that conceals a number of hazards that we wish to 
point out. 



Alternative formulation. Replacing in the original 
Hamiltonian Ti. the spins se with m a or b + 5si and ne- 
glecting terms quadratic in the Ss£S, we obtain our new 
mean-field Hamiltonian 7i^ t : 



TL' m[ = —2N Jm a m,b + AJrrib Sj 



AJnia^Sj 

jeb 



(11) 



In the present formulation, the mean-field approximation 
can be viewed as neglecting correlations between near- 
est neighbour spin fluctuations. The difference between 
the above TL'^ and the 7i mf that appears in (J2J) lies in 
the additional constant term —2NJm a rrib that features 
a temperature dependence through the magnetizations 
m a and nib that must carefully be kept track of. It is 
easily checked that following the above variational pro- 
cedure route 1 with H'^ instead of H mf leads to the same 
results 0. At this stage, another route can be followed, 
that differs from the variational procedure. We decom- 
pose this second route into two steps. 

Route 2a: Self- consistency. From the mean-field 
Hamiltonian 7i' nf it is easy to deduce both the mean-field 
partition function Z'^ and the average magnetization. 
We find 

Z' t = 2 N e 2N(3Jmam " (cosh(4/3 Jm a ) cosh(4 /3 J m h )) N - 12 

(12) 

and 



i a = (s iea ) = jrJ2 s *e~ mLs = -tanh(4 Jf3m b ) (13) 
m b = (sjeb) = -tanh(4J/3m a )(14) 



This system of equations is exactly the one found in llOfl . 
The self-consistency equations (|13f) and (|14f> have only 
the paramagnetic solution when T > Tjv , while a nonzero 
solution continuously develops as T is lowered below Tjv- 
It is then argued that below Tjy, which is identified with 
the Neel temperature, the m a = nib = solution is unsta- 
ble while the solution m a = —nib ^ becomes stable and 
is the physically relevant one. Either more precise discus- 
sions about stability issues are discarded or one finds in 
standard textbooks the following assertion to justify the 
choice of the nonzero solution below Tjv : It becomes sta- 
ble below the Neel temperature (this is true), as can be 
checked by studying the minima of the free energy. This 
is this last sentence that we would now like to discuss. 

Route 2b: Free energy landscape. From the expression 
of the partition function given in 112|) one can easily de- 
duce an expression for the free energy 6' (m a ,nib) as a 
function of the magnetizations on the two sub-lattices: 

6'(m a ,m b ) = - k B TlogZ' ml = N[ - 2Jm a m b 



In [4 cosh(4/3 Jm a ) cosh(4/3 Jm b )} 

2/3 

(15) 

We now express that we search for the states that mini- 
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mize the free energy <j>' 



dm a dmi 



= 



(16) 



Within the framework of the simpler ferromagnetic case, 
this is precisely the wording adopted e.g. in ;2]. It is then 
usually commented that the equations (|16fl are equiva- 
lent to those obtained by resorting directly to the self- 
consistency conditions. A plot of the free energy land- 
scape as a function of the order parameter usually follows. 
And indeed for T < Tjv it may be seen that the paramag- 
netic state becomes a global maximum of </>' considered 
as a function of independent variables (m a ,nib). 

However, below the Neel temperature, the nontrivial 
state (m a , nib) deduced from ifHfl) and (fHE|l is simply nei- 
ther a local nor a global minimum of the free energy 
4>'{m a ,m b ): 

(i) there exist other states, at the boundaries of the mag- 
netization domain, that have a lower free energy; 

(ii) (m a , nib) as given by the nonzero solution of (|13fl and 
HI 41 does not even correspond to a local minimum. 

It is instructive to examine the shape of the free energy 
landscape as a function of the order parameter compo- 
nents (m a ,nib), as plotted in Figure We find that 
the state that globally minimizes the free energy <j>' is 
the fully ordered ferromagnetic m a = nib = 1 state (or 
equivalently m a = mb = —1), whatever < T < Tat. 
Furthermore, the anti-ferromagnetic state corresponds to 
a saddle point of the free energy landscape. This is best 
appreciated on Figure^] With the chosen parameters the 
correct anti-ferromagnetic state has m a ~ 0.77. We have 
clearly come across an unexpected hazard of the mean- 
field approximation. Finally, note that by artificially di- 
viding a regular ferromagnetic Ising model on a square 
lattice into two sub-lattices with independent average 
magnetizations, one would come across the same kind 
of problems for the mean-field free energy (the ferromag- 
netic state would become a saddle-point, the global min- 
imum would correspond to the fully anti-ferromagnetic 
state, etc.). 

Discussion. We now come back to the variational for- 
mulation and we wish to underline that when one evalu- 
ates 4> at its global minimum one finds that 



{m ai m b ) = <j)'(m a ,m b ) 



(17) 



where m a and nib are the functions of temperature solu- 
tion to the system (|13I14|) or JTDJl. This tells us that it is 
perfectly legitimate to follow route 2a using the numeri- 
cal value of <f>' at the values given by I|13I14|) for finding 
the physical solution to the problem. In other words, 
route 2a with the further computation of tj> for the self- 
consistent magnetizations is correct. However it should 



be prohibited to freely vary the magnetizations m a and 
mb as in route 2b to study the free-energy </>' landscape 
and to rely on the latter landscape to discuss stability 
issues. There is indeed no variational principle under- 
lying the derivation of <p' . This is at variance with the 
safe route 1 relying of cj>. This also means that route 2 
can never be used to discuss metastability issues, even 
restricting to self-consistent magnetizations. 

The inconsistency we have brought forth for route 2b 
appears only when the mean-field order parameter is not 
a scalar (as is the case in the ferromagnetic model): Or- 
dering phenomena on more complex substructures of the 
original lattice would inevitably lead to similar results. 
The key point in our second formulation of the mean- 
field approximation (route 2) is that by their very defini- 
tion, m a and nib are the average magnetizations: They 
cannot be considered as freely varying variables. They 
are functions of the temperature determined by the self- 
consistency equations. It is curious to note that the phys- 
ical solution (m a ,mb) to the problem (which is by defi- 
nition a minimum of <p) corresponds to a saddle point of 
the function </>'. 

Therefore we would like to conclude by warning that 
there is in principle no physical meaning to the mean- 
field free energy 4>' seen as a function of a freely varying 
order parameter. Comparing mean-field free energies is 
meaningful only for solutions of the self-consistency equa- 
tions. We have shown this below the Neel temperature, 
but similar problems arise in the high-temperature limit. 
Indeed, expanding the free energy Eq. i|15fl in the vicinity 
of j3 = yields (j)' mf ~ —2JNm a m,b from which one could 
be tempted to conclude erroneously that the stable state 
is fully ferromagnetic when it is of course paramagnetic! 
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FIG. 1: A two dimensional anti-ferromagnet, for T <C Tjv- 
The spins align almost perfectly in opposite directions on two 
sub-lattices a and b of the original lattice. 
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